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Abstract 

We consider possible scale-dependence of the non-linearity parameter /nl in local and quasi- 
local models of non-Gaussian primordial density perturbations. In the simplest model where the 
primordial perturbations are a quadratic local function of a single Gaussian field then /nl is scale- 
independent by construction. However scale-dependence can arise due to either a local function 
of more than one Gaussian field, or due to non-linear evolution of modes after horizon-exit during 
inflation. We show that the scale dependence of /nl is typically first order in slow-roll. For some 
models this may be observable with experiments such as Planck provided that /nl is close to the 
current observational bounds. 



1 Introduction 

Primordial density perturbations are traditionally described by a Gaussian distribution, characterised by 
an almost scale-invariant power spectrum. However the detailed information about the primordial density 
perturbations over a range of cosmological scales offers the opportunity to test in detail the nature of the 
primordial perturbations, both their scale-dependence and Gaussianity [1]. 

The local model for no n- Gaussianity has proved a remarkably popular description of possible deviations 
from a purely Gaussian distribution of primordial perturbations. In the simplest case the primordial New- 
tonian potential includes a contribution from both the local value of a linear Gaussian field and a quadratic 
term proportional to the square of the local value of the Gaussian field: 

= <p G {ft + fhT 1 - (<&)) , (1) 

where /nl^ 1 is a dimensionless parameter characterising the deviations from Gaussianity [2]. (tp^) denotes 
the ensemble average, or equivalently the spatial average in a statistically homogeneous distribution. Note 
that following Komatsu and Spergel [2] we adopt a sign convention for the Newtonian potential $ which is 
the opposite of that used, e.g., by Mukhanov et al [3]. 

In Fourier space we define the power spectrum and bispectrum as 

= (2nfP, s ,(k 1 )S 3 (h + k 2 ) , (2) 
(®k®k 2 *k 3 ) = (27r) 3 J B$(fc 1 , fc 2 , fc 3 )<5 3 (fci + h + k 3 ) , (3) 



1 E-mail: byrnes@physik.uni-bielefeld.de 

2 E-mail: s.nurmi@thphys.uni-heidelberg.de 

3 E-mail: g.tasinato@thphys.uni-heidelberg.de 

4 E-mail: david.wands@port.ac.uk 



1 



and the amplitude of the bispectrum relative to the power spectrum is conventionally given by the non- 
linearity parameter 

f /, , k s = B <s ,(ki,k 2 ,k 3 ) 

X ' 2 ' 3> ~ 2{P*(k 1 )Pz(k 2 ) + Pz(k 2 )Pz(k 3 ) + Pt > (k 3 )P <s> (k 1 )}- 1 > 

In the special case of the local model (1) we have /nl = /nl an( i ^ ^ s clear that /nl is, by construction, a 
constant parameter independent of spatial position or scale. 

This local model turns out to be a very good description of non-Gaussianity in some simple physical 
models for the origin of structure. In particular it can describe the primordial density perturbation on 
super-Hubble scales predicted, up to second-order, using the (5iV-formalism [4-6] if the local integrated 
expansion is a function of a Gaussian random field, N(a g ), at some initial time, <r g (x) = a(ti,x). We then 
have 

S = |C = • (5) 

A good example is provided by the simplest curvaton scenario [7-9], where quantum fluctuations of a 
weakly-coupled field during inflation are well described by a Gaussian random field and the primordial 
density perturbation is determined by the density of the curvaton field when it decays, which is proportional 
the the square of the initial local field [10, 11]. 

In this paper we consider extensions of the simple local model (1). In particular we will characterise 
and quantify the scalc-dcpcndcncc of the parameter /nl which arise in realistic inflationary models for the 
origin of structure. Scale-dependence arises due to two key features. Firstly we examine a multi-variate 
local model where the local expansion is a quadratic function of more than one Gaussian random field. In 
this case scale-dependent /nl can arise if the Gaussian fields have differing scale-dependence, leading to a 
change with scale in the correlation of quadratic terms with the linear perturbation. Secondly, we show that 
scale-dependent /nl can arise even when the expansion is a function of a single canonical scalar field. In 
this case, it is due to the development of intrinsic non-Gaussianity associated with the non-linear evolution 
of the initially Gaussian fluctuations after Hubble-exit. We will refer to this case as a quasi-local model. 

In the case of the so-called equilateral /nl, which arises for example from DBI inflation (see e.g. [12]), the 
scale dependence has already been quite well studied both in theoretical models [13-16] and forecasts have 
been made for future observational prospects [17, 18]. For a discussion of the different possible shape depen- 
dences of the bispectrum see [19]. However for local-type models there has been little previous consideration 
of a scale dependence. The scale dependence of local-type /nl was first calculated in Byrnes et al [20] for a 
specific model of hybrid inflation with two-fields. It was found that although the scale dependence is slow-roll 
suppressed, it depends on a particular combination of slow-roll parameters which is not negligible and can 
easily be larger than the spectral index of the power spectrum. The scale dependence of /nl was considered 
in the case of an ekpyrotic universe in [14]. In the case of an exact solution of two-field inflation, which 
can give rise to a large non-Gaussianity, it was shown that /nl is scale independent [15]. The observational 
prospects for local-type models were considered in [18] which showed that the CMB data is sensitive to a 
scale dependence of /nl- However they used a very simple Ansatz for the scale-dependent /nl: here we 
calculate for the first time the full scale dependence of a scale-dependent quasi- local /nl- 

We note that higher-order contributions to the primordial perturbations are expected beyond quadratic 
order. These only affect the bispectrum at subleading order in the scalar field perturbations, although in 
some cases they might provide the dominant contribution to observables [21]. In the language of [22] they 
are loop corrections. The effective scale-dependence of /nl due to higher-order terms is examined explicitly 
in [23] and was considered previously, for example, in [24]. Notice that the results depend at leading order 
on the infra-red cut-off, and that no infra-red complete theory is yet known [25]. It would be interesting 
to further develop this interesting issue in a case in which the infra-red effects arc well-understood. Note 
that the loop term is not well described as having constant scale dependence, unlike the cases we consider 
in this paper. We also note that [23] only consider the scale dependence coming from the logarithmic term 
which depends on the cut-off, while the other terms which multiply this will also have a scale dependence in 
general. This scale dependence can be calculated using the methods presented in this paper. 

This paper is organized as follow. In Section 2, we focus on multi-variate local models, which are models 
that contain contributions to the primordial curvature perturbation from more than one Gaussian field, in 
which the scale dependence of /nl is due to the different scale dependences of the fields that drive local 
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expansion. As an example of this we study a mixed inflaton and curvaton scenario in sec. 2.1. In Section 3, 
we show that scale dependence of /nl can arise also for single field models, due to the non-linear evolution of 
initially Gaussian fluctuations after Hubble exit. In Section 4 we extend the discussion to systems involving 
multiple fields, taking into account the non-linear evolution of fluctuations after horizon crossing in this 
context. We conclude in sec. 5. 



2 Multi-variate local model 

We will begin by showing how scale dependent /nl can arise in an idealised model that can describe per- 
turbations from multi-field inflation. A more systematic discussion on the scale dependence of this quantity, 
taking into account further second order-effects not considered here, will be developed in later sections. 

A natural extension to the simplest local model given by Eq. (1) comes from considering a local function 
of more than one Gaussian field: 

= E^c® + E fu {VgWVgW ~ (¥>gVg» ■ (6) 
i i, j 

where ip 1 (x) describes n independent Gaussian random fields and fu = fji are n(n + l)/2 constant param- 
eters. 

In Fourier space we have 

(tpltpi) = (27r) 3 Pj(k)S IJ 5 3 (k + k') (7) 

and all higher-order moments vanish for Gaussian fields. It is thus straightforward to construct the power 
spectrum 

P^k)=J2 P vi(k), (8) 
i 

and the bispectrum 

B^(k u k 2 ,k 3 ) = 2^f IJ [P vI (k 1 )P vJ (k 2 ) + P vI (k 2 )P vJ (k 3 ) + PyMPyjih)] . (9) 
i,j 

Thus the non-linearity parameter defined in Eq. (4) is 

Ei.jfij[P V i(ki)P v >j(k2)+P v i(k 2 )P v j(k 3 ) + P vI (k 3 )P v j(k 1 )} 

/NLl*i,*a,*3j [P $ (fc 1 )P$(fe) + P$(fe)P*(fc3) + P*(fc 3 )P$(fci)] ■ { j 

If we restrict our attention to equilateral triangles with ki = k 2 = k 3 = k we have 

i,.j 

where the weight given to the non-linearity of each field depends on their contribution to the total power 
spectrum 

P v i{k) 



-^w- (12) 



The scale dependence is then given by 
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_ dln|/ NL | 

nf - - (13) 

Hlj,k( t i + t J - 2r g )w J (fc)w J (fc)w Jf (fc)/ J j 



J2i,j w i( k ) w j( k )fu 

S/,j( T / + T j)wi{k)wj{k)f u 
/nl 



2(n - 1) , (15) 



5 Notice that we follow [18] but not most previous papers in denning Tif NL = (as opposed to —1) as corresponding to a 
scale independent /nl- However we use the notation Tif rather than njvG since this can be easily generalised to consider a 
scale dependence of the trispectrum. 
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where 77 = d In P v i / d In k and the total spectral tilt is 



dlnP$ „ dlnPr „ v-^ „ . 

n -^^ +3 = ^J +i = ^ WlT > + i - (16) 

Although we have here calculated nf NL in the case of an equilateral triangle, we discuss in sections 3.3 and 
B how to treat the more general case. Note that if all the fields share the same spectral index, 77 = n — 4 
for all /, then nf NL = 0. However we will see in Sec. 4 that /nl could also acquire a scale dependence even 
if all of the 77 are the same, when non- linear evolution of second order fluctuations spoils the Ansatz (6). 



2.1 Mixed inflaton and curvaton scenario 

As an example we consider an idealised model of an inflaton, </>, plus curvaton, a, whose large-scale pertur- 
bations at some fixed initial time, after Hubble-exit during inflation, can be described by Gaussian random 
fields. The mixed inflaton-curvaton scenario has previously been studied by several authors [26-30]. Pri- 
mordial density perturbations due to adiabatic fluctuations in the inflaton field during inflation, ip^, arc 
necessarily very close to Gaussian, = 0, and independent of the curvaton perturbations, = 0. 
However primordial density perturbations due to isocurvature curvaton perturbations may have significant 
non-Gaussianity [11], f aa 7^ 0. We thus have a bi-variate local model 

= <p 4 (x) + <p a (x) + / (7(j (<pl(x) - (cpl)) . (17) 

The resulting non-linearity parameter for equilateral triangles (11) is 

/ NL (fc) = wl{k)f aa , (18) 

where the scale-dependence arises solely due to the scale-dependence of the weighting function, that is given 
by w a (k) = P a (k)/P$(k). We find (here r a = d\nP a /d\nk) 

n fNL = 2(r CT + 3) - 2(n - 1) . (19) 

Using the the results of Wands et al [31] for the scale dependence of the inflaton and curvaton perturbations, 
we obtain 

n-1 = - (6 - Aw a ) e + 2(1 - w a )r\^ + 2w a r\ aa , (20) 

n f NL = 4 (1 - w <?) ( 2e + Vera - V4>4>) (21) 

where, in the notation of [31], r] aa , and e are the usual slow-roll parameters and w a = cos 2 A denotes the 
correlation between the curvaton perturbation, ip a , and the primordial curvature perturbation, $. Note that 
in the limit in which curvaton perturbations dominate over the inflaton in the primordial density perturbation 
(w a = 1) the scale-dependence vanishes. 

It is interesting to estimate the size of nf NL compared to the spectral index. Since both are a function of 
three slow-roll parameters as well as w a , it is not possible to make any definite relation in general. But in 
small field models of inflation one has e <gC {tj^I and it is reasonable to assume as well that r\ aa <C \n<t><t>\- Note 
that since the quadratic curvaton models requires r\ aa > we would otherwise need an unlikely cancellation 
between the two 77 parameters in order to have a red spectral index in agreement with observations. In the 
case that rj^ provides the dominant contribution to both (20) and (21) we have independently of w a 

n fNL = -2(n-l). (22) 

Observations [32] tell us that this implies nf NL ~ 0.1 which interestingly is roughly at the border of de- 
tectability with Planck, for a fiducial value /nl = 50 [18]. It should be possible to obscrvationally test the 
two conditions which are required for (22) to be valid, by considering additional observables. The tensor-to- 
scalar ratio is given by r = 16e(l — w a ) [31], and using the S N formalism it is possible to relate one term 
of the trispectrum with the bispectrum, as t^l = (6/nl/5) 2 /w,t (using the notation and formulas of [33]). 
Hence there are in principle four separate observables which allow one to individually identify the values of 
the four model parameters in (20) and (21). Finally we note that gNL, which parameterises a different term 
of the trispectrum, does not give an observationally competitive signature in this scenario [28] . 



4 



3 Single field quasi-local models 



We now move towards developing a systematic formalism to compute the scale-dependence of /nl m infla- 
tionary models. We start by considering the simplest case where the primordial curvature perturbation is 
generated from fluctuations in one scalar field. We consider two examples. In the first one, we discuss a 
standard single slowly rolling inflaton field with canonical dynamics. While the non-Gaussianities generated 
in this model are unobservably small [34] , this case serves as a useful example illustrating the physical origin 
for the scale dependence of /nl- As a second example, we consider a curvaton scenario where the effect 
of inflaton perturbations can be neglected. This can generate observable non-Gaussianities. We show that 
while the simplest quadratic curvaton scenario predicts a scale- independent /nl ; models with interaction 
terms in the curvaton potential typically give rise to scale-dependence. 
We assume the scalar field has the canonical Lagrangian 



C = - v(</>) 



and that it is light during the inflationary epoch: V" <C H 2 . We assume the field obeys slow-roll dynamics 
during inflation but we do not require its energy density to dominate the universe. We set mpi = (SttG)" 1 ^ 2 = 
1 and introduce the slow-roll parameters 

__A -IfZ!_V -Y!L * _ v'"v m , 

which are all assumed to be small during inflation. If the energy density of the scalar field (f> dominates the 
universe during inflation, we can equate = e, but in a more general set-up these two parameters can differ 
even at leading order in slow roll. 

We will use the SN formalism to characterise the primordial curvature perturbation, given by (5), and 
analyze the scale dependence of /nl- The curvature perturbation on supcrhorizon scales is given by the 
expansion 

C(x) = Nt(ti)fi<Kti,2) + \n^{U) {H(U^f ~ <<ty(*i,z) 2 » + ... , (24) 

where N(U) denotes the number of e-foldings from an initial spatially flat hypersurface at ti to a final uniform 
energy density hypersurface, which we assume is chosen such that the curvature perturbation has frozen to 
a constant value £ = 0. and denote derivatives with respect to the initial value of the scalar field 
4>(U). For the actual computations, it is convenient to transform (24) into Fourier space writing 

C e = N^ti) <%(i,) + ^N^iU) (5<P * {U) + . . . , (25) 

where * denotes convolution and subtraction of the zero mode C^ = o horn (25) is implicitly understood. 
The primordial power spectrum is given to leading order by 

P ( (k) = N$(ti)Ps*{ti,k). (26) 

The calculation of the bispectrum and the derivation of the general form for the /nl parameter, defined in 
eq. (4), are straightforward. One obtains 



-fMh,k 2 ,k 3 ) - "^y+ P c(W (fe) + 2perms • W 



The second term on the right is proportional to the connected part of the scalar field bispectrum -B^, which 
is the three-point correlator of the scalar perturbations. 

The initial time ti in (24) or (25) can be chosen to be any time after the horizon exit of a given mode, 
U > t*(k) where i*(fc) is determined by k = a(t*)H(t*). The curvature perturbation is by construction 
independent of the choice of the initial time ti [35], as discussed in more detail in Appendix A. However, 
statistical properties of the scalar field perturbations 8cf)^(ti) do depend on the choice of U. 
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For a canonical scalar field, slowly rolling during inflation, the intrinsic non-Gaussianity is slow-roll 
suppressed at Hubble exit [34,36]. In our quasi-local model we will thus set the scalar field bispectrum 
B&<t>{ki, k2,ks,ti) — when tj = t*(fci) for all fej. However at later times the distribution in general develops 
non-Gaussian features because of the non-linearities of field equations, and it is this that leads to scale- 
dependence of /nl- 



3.1 Inflaton field 



The case of single field slow-roll inflation is particularly instructive, since there exist general formulas con- 
necting the derivatives of N with slow-roll parameters. Moreover, it is possible to appreciate the connection 
between the scale dependence of /nl and evolution of the second-order curvature perturbations after horizon 
exit. 

Before focusing on the properties of /nlj we recall how in the familiar case of a single inflaton field 
the scale dependence of the power spectrum can be computed in the SN formalism in two different ways. 
Understanding this familiar example will be useful to tackle the more complicated case of the scale dependence 
of /nl- 

The solution of the first order equation of motion for the inflaton perturbations at a time t > t*(k) soon 
after horizon exit can be expressed as [37] 



W) = 



iH(t) 



1 



c + In 



a(t)H(t) 



(3e - r)w) 



(28) 



where e = £0 for the inflaton, is a classical random variable satisfying (a^g) = (2ir) 3 6 3 (k + k'), and 
c = 2 — In 2 — 7, with 7 being the Euler-Mascheroni constant. To leading order in slow roll, the time 
dependence of the slow-roll parameters defined in (23) can be neglected. 

Using equations (25) and (28), one can work out the dimensionless power spectrum of £, defined by 
(Cfej^Cfe ) — (27r) 3 i5(fci + k^li^V c {k\)/k\. To leading order in perturbations, and next-to-leading order in the 
slow-roll expansion, this gives 



V c {k) = Nl{t i ) 



4tt 2 



1 +2£ 



111 



a{U)H{U) 
k 



(3e — 77^) 



and the spectral index 



1 



din 7? 



C 



dlnfc 



can be immediately computed by differentiating (29) with respect to In k. This yields, 

n — 1 = 2?;00 - 6e . 



(29) 



(30) 



(31) 



In practice the same result for the spectral index is usually derived in a different way, that will turn out 
to be useful for the following discussion. In deriving (29) we have not specified the arbitrary initial time U 
in (25), except for the technical constraint that ti needs to be soon after the horizon exit, ti > t*(k), for the 
solution (28) to be valid. A particularly simple choice of ti is to set it equal to the horizon crossing for each 
mode separately U(k) = t*(k). 

To obtain the spectral index we can use the identity [5] 



H 







dh\k 



d_ 

Of, 



k/a 



(32) 



applied to eq. (30), choosing ti at the epoch of Hubble crossing. The second term in the right hand side of 
(32) gives zero when applied to the power spectrum, due to the latter being independent of ti at fixed k. 
Therefore, at leading order in the slow-roll expansion, we have 



I = 2 



dln|Ag 
HdU 

2 ri0 - 2e , 



; dlniJ 
1 HdU 



(33) 
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after defining (and that will be useful in the following) as 

dln|7V | 
"* = ~Hdt~ ' 



(34) 



HdU ■ (35) 

For the inflaton we have N$ = —H/<f> and hence 

n<f> = r] - 2e . (36) 

Thus the expression (33) reduces to the usual formula (31) when focusing on single field slow-roll inflation. 

Similarly, we can calculate /nl in (27) and its scale-dependence using the fact that the curvature pertur- 
bation is independent of ti, sec Appendix A, provided that ti is chosen after each of the scales kj, j — 1, 2, 3 
has crossed its Hubble scale. We concentrate here on equilateral triangles, fc, = k postponing the detailed 
discussion of non-equilateral configurations until Appendix B (but see also Section 3.3). 

If the field perturbations 6<frk are Gaussian at horizon crossing, that is the bispectrum -B£^(fc, k, k, 4*) 
vanishes for any k = a(t*)H(t*), then the second term in (27) vanishes at ti = t*. With this choice, all the 
scale dependence is encoded in the first term of (27), which acquires an implicit dependence on k through 
the function ti = t*(k). Note that for the inflaton we have d/Hdti = (<f>/H)d/d(f> and hence 

n* = ■ (37) 



Thus, using (36) we have 



and the scale dependence is given by 



lhL = ^ = 2e- V . (38) 



din 



n fNL 



- 2n<f, . (39) 



HdU 

To leading order in slow-roll, the result can be written explicitly as 

n fNL = v _ 2e • (40) 

This shows that for single field inflation rif NL is fully determined by the slow-roll parameters evaluated at 
horizon crossing and is first order in slow-roll. (Note that for r\ — 2e = the logarithmic scale-dependence 
diverges simply because /nl is zero.) 

We can physically understand the result by analyzing the evolution of second order perturbations in this 
set-up, also using the results of Appendix A. For single field inflation, we can interpret the scale dependence 
of /nl as being due to the evolution of second-order perturbations right after horizon exit. Indeed, by means 
of the 5N formalism, splitting scalar perturbations into first and second order parts as 5(f) = 6i<f>+ ^S2(f>, we 
can expand the curvature perturbation up to second order in perturbations as 

Cj = Ms + l^f 0W £ ) a + (41) 

The first two terms contain only the Gaussian first order perturbation Si4>. The last term, 52<t>, is the 
second order contribution to the curvature perturbation that, evolving after the epoch of horizon crossing t*, 
is responsible for the scale dependence of local /nl in the single inflaton case. In Appendix A we determine 
the behavior of 62(f), expanding at leading order in slow roll and around i*. Assuming that 82(f) vanishes at 
t*, in accordance with 8(f> being Gaussian at this time, we find 

S 2 (P(t) = (8e 2 - 667, + e) (Si<P*) 2 + ■■■ ■ (42) 

V 2e 
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We see that if <520 docs not evolve at the time of horizon crossing, then nf NL in eq. (40) vanishes as expected. 

Indeed, instead of setting fj = t*(k) in (27) we could have derived the result (40) by treating ti as some 
arbitrary time soon after the horizon exit U > t*(fc). Even though we have assumed the perturbations are 
Gaussian at horizon-crossing, equation (42) explicitly shows that non-linearities in the field evolution result 
in a non- vanishing value for at later times, U > t*. Since the second term in (27), proportional to Bgj, 
depends explicitly on the scale k, this makes /nl scale-dependent. For singe field inflation, this is also the 
only source of scale dependence since the first term in (27) is independent on time and, with the choice 
ii > i*(fc), does not depend on the scale. While this discussion provides a transparent physical explanation 
for the origin of the scale-dependence, in actual computations it is convenient to set ti = t*(fc). These two 
choices are perfectly equivalent, as discussed in Appendix B in more detail. 

It is important to emphasize that in this subsection devoted to adiabatic perturbations of a slowly-rolling 
inflaton field we have neglected non-Gaussianity due to the non-linearity of the inflaton field at horizon-exit 
which is of the same order as that which arises due to non- linear evolution after horizon-exit. In any case 
/nl in this model is suppressed by slow-roll parameters and thus we now turn our attention to models where 
non-adiabatic field perturbations can generate much larger non-Gaussianity after horizon-exit [30]. 



3.2 Curvaton scenario 

In the curvaton scenario, primordial perturbations are generated by a curvaton field a which is light and 
subdominant during inflation, and acquires a nearly scale-invariant spectrum of perturbations. After the 
end of inflation, the curvaton energy density dilutes slower than the dominant radiation component. Its 
contribution to the total energy density therefore increases and its perturbations start to affect the expansion 
history, which generates curvature perturbations. The curvature perturbations evolve in time, f ^ 0, until 
the curvaton eventually decays and the decay products thermalise with the rest of the universe. If the 
curvaton is still subdominant at the time of decay, the perturbations can be highly non-Gaussian. Here we 
concentrate on this limit and assume that the universe remains effectively radiation dominated until the 
curvaton decay, a cx t 1 / 2 . In addition, we assume the primordial perturbation is entirely generated by the 
curvaton fluctuations. In this limit the curvaton scenario effectively reduces to a single field model and it 
can be analysed using the formalism developed above. 

We consider a generic curvaton potential containing a quadratic term plus self-interactions [29,38-42], 

1 rr n+i 

V» = -mV + A— , (43) 

where A is some coupling and n > is an even integer to keep the potential bounded from below. We 
assume the curvaton is long-lived enough so that the energy density stored in the interaction term eventually 
dilutes away and the curvaton ends up oscillating in an effectively quadratic potential before decaying. The 
curvaton energy density at this stage is approximatively given by 

where a{t, cr*) is the envelope of oscillations in the quadratic regime. It is proportional to i~ 3 / 4 and to the 
function cr OS c (cr* ) that depends on the initial value of the curvaton field cr* . While cr osc is a linear function of 
cr* for a purely quadratic curvaton potential, it can acquire a highly non-linear form when self-interactions 
are switched on. As a representative example, at the end of this section we will discuss the curvaton model 
with a quartic self- interaction term in the potential and the corresponding form of the function a osc (a*). 

We assume the curvaton decays instantaneously at H = T, where T is the effective curvaton decay rate, 
and denote its contribution to the total energy density at this time by 

■^osc(^) «1. (45) 



dec 4 



Setting ti = t*(k) in (25) and working to leading order in r, one finds that the derivatives of TV in (25) are 
given by [6] 
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where the primes denote derivatives with respect to ct». Using (46), the non-linearity parameter /nl (27) 
becomes 

5 N„ c _ 5 / go BC <CA a? x 
/NL -6^vJ-47( v 1+ K^FJ ■ (47) 

As before, we have assumed the curvaton perturbations are Gaussian at horizon crossing. 

The scale dependence of /nl can be computed using equation (39) with <fi replaced by a. To find the 
quantities n a and n aa , defined by (34) and (35), it is convenient to first express the above results in terms 
of a reference curvaton value cti(ct*, t*(k)) measured at some time t\ during inflation but after the horizon 
exit of all relevant modes, t\ > t*(k). Since ct sc(cti), we can use the chain rule to write 



chTosc(o'i) _ da osc da 1 _ dcr osc /\ / da\\ 



dai 

di* dai di* dai \ * V ) t ' V ®t* 



(48) 



Moreover, as <i is a time event during inflation, the curvaton evolution until this time is governed by the 
slow-roll dynamics 

r>CTl dcr [ ll dt 



(49) 



from which we get 

'dai \ V'(ai) fdaA V'(ai) 



V'(a) J u 3H(t) ' 

V'trrA 

(50) 



daj u V'(a*) ' \dUj at 3H(U) ' 
Substituting these into (48) and using 3Ha* = —V'(a*), we arrive at the result 

fe-'^T- 1 ' (51) 

to leading order precision in slow-roll. This derivation can be immediately generalized for an arbitrary 
function of the form /(fi), and we arrive at the useful generic result d/(eri)/dlogfc = 0. 

Applying these results to (46), we can easily compute the parameters n a and n aa . According to the 
discussion above, a osc and r do not depend on the wavenumber k. Similarly, using the chain rule, we 
can write a' osc — (dai/ da*)(da osc / da{) where the second term depends only on a\ and is constant under 
differentiation with respect to k. Therefore, we have 

dln< sc d f dcr i\ rr ( . w dinner*) 

= H(U) — In — = -H{t,) = r) aa , (52) 



dlog k di* \ da* J dt 

where we have used (50) and the slow roll equation 3Ha = —V 1 . Similarly, we find 



Substituting these into (39), we find the scale dependence of /nl is given by 



V"'(a*) ( er osc er osc 



osc 



It can be immediately seen that nf NL = for a purely quadratic curvaton model. However, in the presence 
of interactions, V'" ^= 0, the non-linearity parameter becomes scale-dependent. Equation (54) gives a rough 
estimate |rif NL | ~ \r}oo\ but the precise prediction depends on the curvaton potential and the details of the 
dynamics encoded in the derivatives a' osc and a" sc . 

As a specific example, we estimate rif NL for the curvaton potential V = 1/2 m 2 a 2 + Act 4 . This case can 
be treated analytically using the results of [42] . In the limit where the interaction term dominates initially, 
Act* > m 2 , the function ct osc (o'*) in (44) can be estimated by equation (4.9) of [42], 

I 3 e -0.80\/i 

aosc(a * } " "* |r(o.75 + zo.5Tyi)| " g * c(s) ' (55) 
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where s = Xal/m 2 . As discussed in [42], this agrees well with numerical results for s > 1 but becomes 
inaccurate when approaching the quadratic limit s = 0. The qualitative behaviour however remains correct 
even in this limit and we can therefore use (55) for arbitrary s, provided that we understand the results in 
the limit s < 1 as order of magnitude estimates only. 
Substituting (55) in (54), we find 



"1k 



£(s) 2 + 2s£(s)£'(s) 



HI VlO^(s)f( S ) + 4 S 2(£( S )£"( S ) + e( S ) 2 ) + £(s) 5 



m 



m = 



24s 
1 + 12s 



(56) 



The function /(s) can be evaluated using (55) and its behaviour is illustrated in Fig. 1. ^From (56) and Fig. 




0.01 



0.1 
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Figure 1: An estimate for the function /(s) in equation (56) shown as a log-linear plot. For s < 1 only a 
qualitative agreement (error ~ 20%) with the exact behaviour is expected. 

1 we see that rif NL vanishes in the quadratic limit s = as expected. If the interaction is initially comparable 
or dominates over the quadratic term in the potential, s > 1, we find rif NL ~ ij aa up to a factor of few. 

Note that for the potential V = l/2m 2 a 2 + Act 4 we have r\ aa > and since n — 1 = — 2e + 2r) crcr , it is 
assumed in [42] that e gives the dominant contribution to the spectral index. In this case it is clear from 
(56) that nf NL will be subdominant compared to the spectral index and hence unobservably small. However 
in the more general case (54) it is possible to have r\ aa < at the time of horizon exit but with a minimum, 
e.g. if the curvaton has an axionic type of potential [43]. In this case it might be possible to have a large 
rif NL . Another case which could lead to a large rif NL in the curvaton scenario is the non-perturbativc curvaton 
decay [44, 45] but an investigation of these points is beyond the scope of this work. 



3.3 Scale-dependence for generic single-field case 

It is straightforward to generalize the result derived for nf NL to any quasi-local model where the curvature 
perturbation is generated from a single field. Indeed, since equation (27) for /nl is valid for all such models, 
we can simply set ti = t* in it and read off the result from equations (34), (35) and (39) where <f> is 
now understood as the generic scalar degree of freedom responsible for generating the primordial curvature 
perturbation. We show in Appendix B that the result rif NL = — 2n^ actually holds not only for the 
equilateral case but for arbitrary triangles that change their area while preserving their shape. 

In Appendix B we also derive a result that allows to study variations that change the shape of the 
triangles. Considering a triangle with sides fci,fc2,&3 not too much different from each other and from a 
generic reference scale fc p , we can expand /nl around k v as 



2 perms 



/NL (fcl ,k 2 ,k 3 ) = / N L 0t y V {n-i, sn-i • ( 57 ) 

Here f^ ot corresponds to /nl evaluated on an equilateral triangle with sides k p and the result is given to 
leading order precision in 0(n — 1, e^, ry^, n^, n^) ln(fcj/fc p ). Note that if is an isocurvature field during 
inflation, like in the curvaton scenario (or more generally in models characterized by a large /nl [30]), then 
e ~Oin (57). 
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The result (57) allows one to generally study the leading order contributions to the scale dependence of 
local /nl^ ^f NL obtained from it reduces to equation (39) for arbitrary triangles that change area preserving 
their shape, while variations that allow the shape to change in general lead to a different scale dependence. 
For example, in the squeezed limit k^ <C k\ ~ k% = k we find the result 

fj \TJ ■ (58) 

If we take a derivative of /nl with respect to the scale k, while keeping the ratio k^/k fixed, we again have 
n fNL = n 4>4> ~ % n d>- But it is clear from (58) that if we would instead take the derivative w.r.t. fc 3 while 
keeping k fixed or vice versa we would get a different result, — — 2e^ + rj^ and — + 2e^ — 
respectively. 

We note that (57) has a more complicated form than the simple Ansatz for a scale dependent /nl used 
by Sefusatti et al [18] (Eq. II. 6 of their paper). The two forms agree in the case of an equilateral triangle 
but in general it is not possible to factorise out the local scale dependence which they denote F(ki, fe, fe). 
It would be interesting to see how the forecasted constraints on nf NL vary if one uses the full form for /nl 
given in (57). 



4 Multi field quasi-local models 



A set-up in which more than one scalar field determines the properties of the curvature perturbation can 
be treated similarly to single field, but in this case there are generically two sources which contribute to a 
scale dependence of the quasi-local /nl- The first, discussed in Section 2, is associated with the fact that the 
various Gaussian fields constituting the system can have different scale dependence, leading to a change with 
scale in the correlation of quadratic terms with the linear perturbations. The second, discussed in Section 
3, is due to the non-linear evolution of second order perturbations after horizon exit. 
The Lagrangian governing the system reads 



C 



1 



Sud^d^ -v(<t> K ). 



(59) 



Notice that for simplicity we assume a flat metric in field space. In this situation, the expansion for the 
curvature perturbation, valid at super-horizon scales, reads in Fourier space 



Cjf(t/) = NMu 



(60) 



where we have truncated the expansion at second order. The quantities Ni, Nu indicate derivatives of the 
number of e-foldings N = J. f Hdt with respect to the scalar fields. In the framework of multi-field inflation, 
they explicitly depend on t/, while the dependence on U is again included in a dependence on <fi L , writing 
Nj = Ni(tf, 4> L {ti)). Again, £(tf,x), as defined in (60), does not depend on tj. 

In order to proceed, it is convenient to know some properties of the solutions for the first order scalar 
perturbations soon after horizon exit, that reads 



(l-e)5 I J + 



In 



a(ti)H(U) 



4 



(61) 



where c = 2 — In 2 
parameters [37] 



7 with 7 the Euler-Mascheroni constant, while e IJ is a combination of slow-roll 

ViVj Vu 



9H 4 'ZH 2 



(62) 



which is small at the epoch of horizon crossing. Finally, a~ is a classical, time independent random variable 
satisfying the orthogonality property (a~ ai 



(2tt) 3 6 ij 5 3 {k + l). 

By proceeding as we did in Section 3, we can calculate the spectral index. By applying identity (32) to 
the power spectrum, one finds 
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n ' 1 = ' 2£+ EjNjNj (63) 
where we introduce the notation (as a natural extension of (34) and (35)) 

dlnjV 7 dlniV/j 

ni = ^m ' niJ = ^idt-> (64) 

and all quantities are evaluated at time U = t* corresponding to horizon exit. We checked that this expression 
coincides with the one of [5], in the case in which the derivatives of N are explicitly expressed in terms of 
slow- roll parameters [46]. 

We can use the above equations and definitions to find an expression for the bispectrum associated with 
the curvature perturbations, 



B Q = (2tt) 3 N r NjN L BI^S^, <ty£, S<f^) 
NiNjNlm (S il + g IJ ) (S JM + g JM ) 

+ N R N S N T N Z (6 TZ + q TZ ) {5 RS + q RS ) 



( P c (fc 1 )P c (fc 2 ) + perms) , (65) 



where the scalar perturbations are evaluated at a common time t j , and Bg^ is proportional to the connected 
part of the three point function for the scalar perturbations. The parameter /nl is obtainable from the 
definition (4), or equivalently 

5 /NL = [P ( {ki)P c (kt) + perms]" (66) 

The tensor q IJ is easily obtained from formula (61). For our purposes, it is enough to say that it can be 
expressed as 

q " = q { J + q i J ln( a -^^\ (67) 



where q± and q2 are linear functions of slow-roll parameters, with scale-independent coefficients. At leading 
order in slow- roll, this quantity depends on k only by means of the logarithm, log k/(aH). 

The quasi-local form of /nl given in eq. (66) has various sources of scale dependence. In order to 
investigate them, we proceed as in the single-field case. We consider scale variations that preserve the 
triangle shape: here we focus on equilateral triangles, while the case of triangles with different shape is 
discussed in Appendix B. In order to calculate nf NL , we apply formula (32) to equation (66). Because /nl 
is independent of ti, the second term of eq (32) does not give contributions. The first term acts only on the 
first term of (66), since when applied on the second term it gives only contributions subdominant in slow-roll 
that can be neglected, analogously to what happens in single field case. Moreover we can neglect its action 
on q IJ , since this provides contributions that are higher order in slow- roll parameters. 

Then, to conclude, the expression for rif NL in this context is 



nr„, = in 

Hdt, 



N I N J N IJ . 
{N L N L f 



(68) 



which is a natural generalization of what happens in the single field case. We make a comparison with 
the explicit multi-variate local model discussed in Section 2, by expanding equation (68) by means of the 
quantities ri\ and n/j defined in eq. (64). One finds 

E/,j("J + "J + n IJ )N I NjN IJ _ 4 Y^niNim 
nfNL El,mN l N m N lm EjNjNj ' 

Ej j(m + nj + nu)NiNjNu 



Elm n lNmN lm 



- 2 (n c - 1 + 2e) . (69) 



Eq. (69) can be seen as the analogue of equation (15), this time calculated by using the SN formalism. If 
not all of the rij are the same, there is a source of scale dependence for /nl, mainly due to the fact that 
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first order, Gaussian perturbations have different scale dependence (see Section 2). But, even if all the fields 
have the same n/, implying nj — ^ (n^ — 1 + 2e), rtf NL still will not generically vanish due to the evolution 
of second-order perturbations after horizon exit (see Section 3). 



4.1 Multi curvaton scenario 

In this section we consider the possibility that two curvaton fields contributed to the primordial curvature 
perturbation. Following [47] we assume that the curvatons are non-interacting and we neglect the inflaton 
field fluctuations (for a related scenario see [48]). In the limit that the non-Gaussianity is large, which is 
the observationally interesting case, we can write the primordial curvature perturbation in the multivariate 
local form 



C = r a ( a + r b Cb - r a ( CGa + 7 CGa + r M CGfc + 7 Cg6 ) . (70) 



where £ a is the curvature perturbation of the curvaton field a and (g a is the Gaussian part of its perturbation. 
We denote the initial, horizon crossing curvaton field values a* and 6* . Then the Gaussian part of the of the 
fluctuation due to the curvaton field a is given by 

t - 2Sa * (t\\ 
CGa - o ( 71 ) 



and it is useful to further define 



Then 



3 a* 



P = T L , r = -- (72) 



K n 



P( = r 2 a P Ca +r 2 b P Cb = (rl + r 2 b p 2 )P Ca , (73) 



4 (r 2 a + P 2 r 2 f 



/nl = 7 ^2 , m r 2^2 ■ (74) 



The full result for /nl is given in [47], but the simplified expression above is a good approximation for 
|/nl| > 1. 

Since this form of the curvature perturbation satisfies Ansatz (6) we could have analysed this model 
using the results of section 2, but we use this model to provide an example of using the formalism of the 
previous section. To calculate the scale dependence we can use (69), and from (70) and (71) we can read off 
the derivatives of TV as 

N a = \-, N aa = 2^, (75) 
3 a* 



and similarly for the derivatives with respect to the b field. Thus, we can proceed as in Section (3.2) to find 
ft a = Vaa, n aa — 2rj aa and similarly for n^. Putting all of this together we find 



9 , 2r ' 2 ^ 2 ^ 



ftfNL = ^(Vaa ~ Vbb) 



^3 



P 2 /? 4 



r 3 +(3 i r 2 + 1 r 3 + /? 4 r 2 + /3 2 



(77) 



Notice that in the limit 2 — > or f3 2 — > oo only one of the fields contributes to £, since P^ b = (3 2 P( a • In 
that case we can see from (77) that rtf NL — > in agreement with the expectation from section 3.2, because 
we are considering the non- interacting case. Also in agreement with Sec. 2 we see that nf NL = if r\ aa = r\bb- 
Since both of the terms of n — 1 which depend on the rj parameters must give a positive contribution to the 
spectral index and we observe a red spectral index the contribution from these terms must be small barring 
unlikely cancellations between these two terms and the — 2e. In addition each of the four fractions in (77) is 
less than unity, so we conclude that (77) is likely to be suppressed compared to the spectral index. 
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5 Conclusion 



We have studied generalisations of the local model of non-Gaussianity, thereby seeing in which cases the 
standard assumption that it can be described by a single, scale-independent parameter is valid. We have 
shown that this is only strictly valid in specific models, where the primordial curvature perturbation is 
generated by a single scalar field which acts as an isocurvature ("test") field during inflation and has a 
quadratic potential. An example of this is the simplest curvaton scenario. However as soon as the scalar 
field has interactions these generate non-Gaussianity of the field perturbations after Hubble exit, and these 
give rise to a scale dependent non-Gaussianity which we call quasi-local, meaning that in the limit that the 
scale dependence of /nl goes to zero one recovers the local model. 

We have also shown that even in a multi-variate local model where more than one Gaussian field con- 
tributes to the primordial curvature perturbation the effective /nl has a scale dependence unless all fields 
have the same scale dependence. An example of this is a mixed inflaton and curvaton scenario or a multi- 
curvaton model, where the curvatons have quadratic potentials with different masses. 

We have developed a formalism, based on the 6N approach, which allows us to obtain compact expressions 
for nf NL = din \ f^\/dlnk, in both the single and multi- field case. We have also discussed more generally the 
validity of defining a scale dependence of /nl with respect to a single scale when the bispectrum is generally 
a function of three variables, except in the case of an equilateral triangle. We have shown that the scale 
dependence of /nl is independent of the shape of the triangle it describes, provided that we consider scale 
variations which preserve its shape. It is then appropriate to perform the simpler calculation of /nl for an 
equilateral triangle, and use this to calculate rif NL by taking its logarithmic scale dependence. 

We applied our formalism to discuss several specific models. Our results suggest that the scale dependence 
is typically first order in slow roll. The precise value however depends quite sensitively on details of the model 
which, makes it in principle possible to use rif NL to discriminate observationally between different models. 

We have shown that while the simplest realisation of the curvaton model has an exactly scale independent 
/nl, almost any extension to a more realistic scenario does give rise to a scale dependence. In the case of an 
interacting curvaton or a multiple-curvaton scenario the scale dependence is likely to be suppressed compared 
to the spectral index. This is because nf NL in these scenarios depends schematically on ?7curvaton which is 
normally positive and, in light of the observational preference for a red spectrum, this is likely to be small. A 
detection of both /nl and rif NL would therefore be a signal of non-trivial dynamics in the curvaton scenario, 
or that the inflaton perturbations have also made a significant contribution to the observed power spectrum. 

In the case of a mixed scenario, in which the primordial curvature perturbation has contributions from 
both the inflaton and curvaton perturbations it is quite natural for there to be a consistency relation between 
the scale dependence of the power spectrum and the bispectrum, nf NL ~ — 2(n — 1), and this could be 
observable with Planck assuming a large enough fiducial value for /nl, close to the current observational 
bounds. In this case the bispectrum would be larger on small scales and this could make a detection of /nl 
using large scale structure data more likely. In any scenario where one includes the Gaussian perturbations 
of the inflaton field as well as the perturbations of a second field which generates non-Gaussianities one will 
have a scale dependent /nl, unless both fields have exactly the same spectral index. 

For the case where non-Gaussianity is generated during slow-roll hybrid inflation one generally has a 
non-negligble rif NL [20]. In this case, since the magnitude of /nl can only grow during inflation, larger scales, 
which exit the horizon earlier, will be more non-Gaussian and one necessarily has rif NL < 0. 

It should be straightforward to extend the formalism we have presented here to the primordial trispectrum 
(the four-point function) [33,49] and we expect that the two non-linearity parameters which model it, t^l and 
9nl, would also have a first order in slow roll scale dependence. We have chosen the notation nf NL in a way 
that is extendible in an obvious manner to study this. For example in the case that a single field generates 
the primordial curvature perturbation one has a non-Gaussianity consistency relation, tjvl = (6/nl/5) 2 . 
Then assuming our formalism can be extended in the obvious way it follows that n TNL = 2nf NL . 
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Appendices 

A Explicit proof that / NL is independent of the initial time 

In this section we show explicitly that /nl is independent of ti. In the first part, we focus on single field 
slow-roll inflation, in which the unique scalar plays the role of inflaton field, and we study the equations of 
motion for the perturbations. In the second part, we prove the same results using general properties of our 
definition of the curvature perturbation. 

In the single field case, ( is conserved on super horizon scales so we do not need to consider a further 
dependence on tf . Then the SN formalism provides the following expression 

C = 7B H+ ^ S2 *- <78) 

All quantities on the right hand side are evaluated at a given time U. It is useful to split the scalar field 
perturbation into first and second order perturbations as 5(f) = Si(j> + \&24>- If we choose ti to coincide with 
Hubble crossing, ti = f*, we have 

|fNL = *?-2e| t , + |f$. (79) 

We are following the notation of [46] and fj?) has a complicated k dependence but it is zero if we take 
^20* = at horizon crossing, a condition which we will assume from now on. We now consider what 
happens if we choose to use a time ti > t* . We first note that 

d(r)-2e) _ „ 2 



dN 



= 6e?y - 8e z - ^ , (80) 



where £ 2 = V V Y is a second order slow-roll parameter. We can get the scalar field time dependence 
from [50]. From (2.13) of that paper it is easy to see that 

= M. [1 + H(ti - U)(2e - rj)] , (81) 

We assume that H (ti — £*), corresponding to the number of e- foldings N between and ti, is much smaller 
than the inverse of the slow-roll parameters. Since we only need to consider about ten e-foldings for obser- 
vational purposes, this is an adequate approximation. 

The evolution of second order perturbations is more complicated, and we need to consider (2.15) of [50], 
at leading order in a slow-roll expansion. We find 

-| (M)" + (V - ^4> = yj (8e 2 - 6e V + £ 2 ) S^ 2 . (82) 

We note that the 8e 2 term above is not present in (2.19) of [50], which is because they use a non-standard 
definition of e. For our purposes it is important to include it. Since we only need to solve this equation at 
leading order in slow-roll, we can neglect the additional time evolution of the slow-roll parameters, as well 
as the evolution of 5i<f> 2 . This makes the integration simple, and we find 

= + 4= (8e 2 - 6er? + £ 2 ) S^ 2 , (83) 
V 2e 



For the reason explained above, it does not matter here whether the terms which multiply N are evaluated 
at ti or f*. Using (66), it is then clear that the time dependence of the slow- roll parameters in (79) as given 
in (80) is compensated by the time dependence of as given by (83), rendering /nl independent of ti. 
The previous discussion can be extended to the case of multi-fields. 
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One can arrive to the same result in a different way, that is easily gcncralizablc the case of multiple scalar 
fields. Recall that the curvature perturbation £ is defined by means of the spatial part of the metric 

9ij = a 2 (t)e 2 ^ 7ij (t,x), (84) 

where a(t) and jij(t,x) indicate background quantities. In [35], it has been rigorously proved that, at 
all orders in the expansion of the perturbations, the curvature perturbation at supcrhorizon scales can be 
expressed in terms of a variation in the number of e-foldings: 

C = M{t f ,t h x)-N{tf,U) = 5N(tf,U,x) (85) 

with 



tf i ftf 



N(t f ,U) = I Hdt, JV(t f ,U,x) = I dt[-?—) . (86) 

Jti Jtt \P + V J U 

In the previous equations, the energy momentum tensor for the system is assumed to correspond to an ideal 
fluid: this condition is satisfied for a multiple field set-up. 

It is simple to see that the curvature perturbation £ given in eq (85) is independent of ti, as long as U 
labels a flat hypersurface. Indeed, one has 

±^ = -H(ti) + \ (-^-)(ti,x). (87) 
oti 3 \p + pj 

On the other hand, as discussed in eqs (18) and following of [35], as long as one is focusing on flat hypersur- 
faces, the following energy conservation equation holds 

a p 



a 3(p + p) 

Using the previous relation, (87) ensures that £ is independent of U, at all orders in the expansion, as we 
wanted to prove. Because the curvature perturbation is independent of U, /nl is also independent of this 
quantity. 

The information that ( is independent of ti, allows also to straightforwardly obtain the equations of motion 
for first and second order scalar perturbations 8\4> and Expanding £ at first order in perturbations, one 
has 

C(x) = JV*(ti)M(*i,2) (89) 
The fact that the right hand side is independent of ti provides the equation (we use the notation of Section 
3) 

(S 1( f>y + n^HS 1( j) = (90) 

leading to cq. (81) for single field slow roll inflation. Analogously, expanding ( up to second order and 
proceeding in the same way, one finds the following equation of motion for second order perturbations, valid 
at leading order in slow roll 

{5 2 <j))' + 7i0 H 5 2 <p = HN^n^ [n^-2n^) (S^) 2 , (91) 

that corresponds to equation (82). It should be possible to extend the present method to the multi-field 
case, even with non-canonical kinetic terms, possibly reproducing the results of [51]. 



B Non-equilateral triangles 

In this appendix, we explain how to calculate the running of /nl for generic triangle configurations. We 
directly discuss a multi-field system, the single field case is easily obtained as special limit of our compu- 
tations. To perform the calculation, it is convenient to set the initial time ti, in the expansion (60) for £k, 
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equal to the horizon crossing time tj = t*(fc). With this choice and using Wick's theorem, the bispectrum 
can be written as 



B c (k 1 ,k 2 ,k 3 ) = A f ( 1) ;JV(2)jiV(3) tM 



dq 

(27r) J 



(^(iO^qfe))^^)^^^)) +2perms , (92) 



since at horizon crossing we set the connected three and four point functions of 8<j) to zero, thereby assuming 
the perturbations are Gaussian at this time. In the previous formula, we introduced the notation N(j\j = 
Ni(tf,tj) with tj corresponding to the horizon crossing epoch for the mode kj. Notice that the two point 
functions in (92) involve arguments evaluated at different times. In general, to first order in slow-roll, from 
(61) it follows that 



Here we have defined a quantity 



M 



hi 



a(tg)H(t a ) 
a(t b )H(t b ) 

VlV m V lm 



e m 



(93) 



(94) 



9H 4 3H 2 ' 

which differs by a factor 2e8LM from clm defined in (62) since we have explicitly included the time evolution 
of H in (93). Using the relation between time and scale at horizon crossing, k = aH, to leading order in 
slow-roll the two-point functions with unequal time arguments can be written as 



M 



M 



In 



1 ** 

fei 



(2^) 3 <5(k! 



e L M 



KlL 



(HUtiW^ih)) (95) 
(6™ + 2ce IM )P i {k l ) 



(27r) 3 ^ 3 )(k! 



2c + ln{ 



±J N {1)R N {1)s (S^ + 2ce^) 



N il)R N {1)s (S^ + 2ce^) 



and similarly for the other terms in (92). Recall that c = 2 — In 2 — 7, with 7 the Euler-Maschcroni constant. 
The result is independent of the time at which the slow-roll matrix i IJ is evaluated since we work to leading 
order in slow- roll. Using (95), the bispectrum (92) can be written as 



B ( (ki, k 2 ,k 3 



N(i)iN(2)jN( 3 ) LM 



2c + ln£^ ) e IL 



JM 



( 2c + lnf|)e 



Nn )R N il)s N i2)T N (2)z (5*s + 2ce««) ($tz + 2c e TZ ) 



'(l)K- iv (l)5- iv (2)T- iv (2) 

xP c (k-L)P c (k 2 ) + 2 perms , 



(96) 



corresponding to equation (65) rewritten using the choice t% = t*(fcj) in (60) for each of the three modes ki. 
From (66) one finds the corresponding non-linearity parameter /nl- 

To evaluate the scale dependence of /nl, we perform an expansion around a pivot scale k p not too 
different from fcj's. Let us start by considering scale variations that preserve the (arbitrary) shape of the 
triangles: we then write the wavenumbers ki in terms of dimensionless parameters on and a common scale 
ki = koti. Here we treat the parameters on as constants, concentrating on variations of the scale k that 
preserve the shape of the triangle. Denoting the derivatives of N at the pivot scale by 



din TV/ 



nu 



d In Njj 



and defining 



we can expand /nl around k p as 



hh(k) = hh(k p ) 



1 + 



Hdt a=k p /H Hdt 
" EjNjNj ' 

'Ej,j(»/ + n J + nu)NiNjNij 
Elm n lN m N lm 



a— k p /H 



4n In 



(97) 



(98) 



(99) 
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Here the coefficient of the logarithm is given to leading order in slow-roll parameters and to the same precision 
the non-linearity parameter at the pivot scale reads /nl(^p) = (^/6)NiNjNjj / (NlNl) 2 . From equation 
(99) we see that for a generic non-equilateral configuration with a fixed shape, the scale dependence of the 
non-linearity parameter is given by 

Y Jh j{ni + nj + n I j)N I NjN IJ 

"fNL = N N N 4n ' (100) 

to leading order precision in slow-roll and in the expansion in In (k/k p ). To this precision, the result coincides 
exactly with the one obtained for the equilateral case, equation (69), and is independent of the constants oii 
parameterising the shape of the triangle. 

However, if both the scale and the shape of the triangle vary simultaneously, the result will in general 
be different. In principle, it is straightforward to work out the scale dependence for an arbitrary variation, 
starting from the general equation (96). As an example we discuss the single field case. The bispectrum (96) 
depends on e/j which in the single field case has only one component e<^ = 2e^ — 77^. Equation (96) can 
be expanded around the pivot scale k p as 



B(k 1 ,k 2 ,k 3 ) = ^ 



1 + (n<f,<f, + 2e^) In ( ^) - (n^ + e^) In (j^-j - (n^ + e^) In ( ^ 



P c (fci)P c (fe) , 



plus permutations. To the same order of accuracy as in (99), we can then expand /nl around the pivot scale 
as 

(H (¥) +2perms 
/ NL (fci, k 2 , k 3 ) = fZ 0t ^ ; ;{ c - 4/ x » c - 4 V W • (101) 

(fe) (fe) + 2 P™ 

This expression allows to study the scale dependence of /nl in full generality for arbitrary variations of the 
triangle. 
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